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ABSTRACT 


The issue of orbital debris resulting from human activities in space is a growing concern for the space 
users' community. Waste generated in Low Earth Orbit (LEO) can stay in orbit for a long time, 
creating significant hazards for other spacecraft flying at lower intercepting orbits. Many spacecraft, 
especially crewed vehicles, are required to vent fluids into space. These fluids include propellant, 
wastewater, excess condensate, and others. It is important to analyze the behavior of particles that 
result from these activities, since each individual particle is capable of damaging or destroying a 
spacecraft in a lower, crossing orbit, and such particles are invisible to tracking radar systems on the 
ground. The deorbit trajectory of an ice particle is complex. It depends on factors including attitude of 
the vehicle during vent, initial velocities of particles, altitude at which the vent occurred, and 
numerous evaporation and sublimation factors. These include contamination within the vented water, 
evolution of bubbles within the clear water, and sun flux factors such as time of the year and current 
beta angle. 

The purpose of this study is to examine the influences of these factors on the trajectories of ice 
particles resulting from condensate water dumps, and to bound the safe trajectories of spacecraft that 
lie below the venting spacecraft. 

NOMENCLATURE 




(in alphabetical order) 


A 

surface area of the particle 

^ therm 

thennalization coefficient 

Alt 

altitude of the particle 

K 

empirical constant 

b n 

ballistic number 

k = 

Boltzmann constant 

c D 

coefficient of drag 

K 

sublimation rate constant 

C p 

ice heat capacity 

L 

heat of sublimation 

C 0’ C g’ C 0 ~ 

surfactant concentrations 

M 

molecular weight of diffusing 
substance 

D 

diffusion coefficient 

m = 

water molecule 
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m c 

= 

complex index of refraction 

N 

= 

GMT day number 

n, ik 

= 

real and imaginary indexes of 

P(T) 

= 

refraction respectively 
vapor pressure above the particle 

P b 

= 

internal pressure of the bubble 

Ps 

= 

pressure at the particle surface 

pM 

= 

internal vapor pressure at T 0 

Q Albedo 

= 

reflected heat flux by Earth 

Q atm. heat 

= 

collisional heat flux 

QEarth 

= 

heat flux from the Earth 

Qin 

= 

incoming heat flux 

Space 

= 

heat flux emitted to space 

Qsunaverage 

average heat flux from the Sun 

Qsunreal 

= 

real heat flux from the Sun 

R 

= 

gas constant 

R„ 

= 

radius of the bubble 

R{) ’ Rdrop 

= 

radius of the particle 

^Earth 

= 

Earth radius 

T 

drop 

= 

particle temperature 

T 

Space 

= 

space background temperature 

V 

= 

velocity of the particle 

dR 

= 

rate of sublimation 

dr 

a{T) 

= 

molecular sticking fraction 

x(R*o„, 


= absorption coefficient 

s 

= 

emissivity coefficient 

Earth 

= 

Earth viewfactor 

P Space 

= 

Space viewfactor 

A 

= 

wavelength 

P atm 

= 

atmospheric density 

P drop 

= 

ice density 

<7 

= 

surface tension 


= 

Stefan - Boltzmann constant 

T 

= 

time 


INTRODUCTION 

When a vented stream starts disintegrating due to 
rapid expansion of dissolved gasses, it creates 
many millions of particles of different size and 
shape. The outer layers of such particles are 
cooled by fast evaporation from the surface. 
Larger particles tend to continue to break up until 
surface tension equilibrates expansion forces that 
otherwise tear particles apart. Previous studies 


predict particle sizes to be comparable to the 
initial diameter of the stream [1] which is up to 
twice the size of the venting nozzle. For this study 
we detennined that droplets will completely 
freeze within 10 seconds after the vent with 
average diameter of 2 mm. Although there is a 
large number of misted particles in addition to the 
main stream, this particular study cannot be used 
for the particles of very small diameters, where 
Mie scattering is applied. We did not perform 
sublimation analysis of such small particles 
because the intent of this study was to characterize 
only the particles that posed a ballistic threat. 
Analysis starts at the time when a droplet is 
completely frozen and continues until the droplet 
reaches an altitude of 250 kilometers, or 
sublimates to one hundredth of its original size. 
Once fonned, droplets cool very quickly and 
reach equilibrium temperature of - 1 90K in the 
order of seconds. This study looks at different 
aspects that affect (A) the sublimation rate, such 
as: beta angles, time of the year, initial altitude, 
impurities, surfactants, etc. and (B) the deorbit 
rate, such as: the initial altitude at which vent has 
occurred, numerous atmospheric parameters and 
the drop/particle cross section, which is 
completely determined by (A) above. 

ASSUMPTIONS 

In order to simplify our study, we made several 
assumptions to our model. We will address the 
rationale and the effect of the assumptions later in 
the paper. For the modeling purposes we assumed 
a spherical crystalline clear ice particle that has no 
impurities or organic solvents. Particles were 
assumed to have a clear surface (no roughening) 
and to be homogeneous (have no dissolved 
bubbles). The background temperature for Earth 
was taken to be 255K at 50 km reference altitude 
and background temperature of space was 
assumed to be 3K. The absorptivity of Earth IR by 
the ice particle was calculated from complex 
indices of refraction to be approximately 
x{R t X) — 0.93 . In the atmospheric model we 
assumed that after venting, the particles remain in 
circular or nearly circular orbit with initial venting 
altitude of 400 km. The coefficient of drag 
C D was assumed to be 2, characteristic of very 
smooth spherical particles. 

MODELING 

In vacuum large diameter streams with high vapor 
pressures have a tendency to burst, forming a 
cloud of particles which subsequently freeze and 
eventually sublimate. The behavior of the stream 
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and the fonnation of the particles in vacuum 
depend on several factors such as vapor pressure, 
viscosity of the vented fluid, dissolved gas 
content, surface tension, surfactants and the 
density. 


smaller particles. Thus there is a lower potential 
for particles to contain bubbles as initial drop size 
decreases. Any pre-existing bubbles will act as 
fracture points in the original stream, not tending 
to remain encapsulated in the liquid. 


PARTICLE FORMATION 

“Boiling” of dissolved gasses and surface heat 
transfer due to rapid evaporation are governed by 
a diffusion equation [2], It can be assumed by 
basic transport considerations that rapid changes 
in a liquid’s temperature or dissolved gas content 
occur in a relatively thin layer near the surface, 
and it will freeze from the outside in, trapping 
most dissolved gasses inside the drop. However, 
we do not expect gasses to evolve into bubbles. It 
is obvious that for a bubble to exist or to grow its 
internal pressure has to be at least greater than the 
competing pressures of the drop and bubble wall 
curvatures: 


cr 2cr 

Pb = Ps + + ■ 


R~drop 


R h 


(l) 


Assuming small thermal gradients in the small 
particle, the internal bubble pressure must be at 
the vapor pressure p v (T 0 ) and p . ~ 0 for a very 
cold surface. 


Thus, a minimum bubble radius that permits 
growth is 


R h =2a(p v --^—)~\ (2) 

R drop 

for such a cold particle surface. Figure 1 
schematically shows bubble growth in the 
particle. 



Figure 1 Bubble expansion in the particle 

From equation (2) we can calculate that for 
cavitations to occur in 2 mm diameter particles the 
radius of the bubble R b has to be greater than 
10% of the radius of the particle, or larger for 


As a particle freezes the expansion of water will 
apply increasing internal pressure to the inner 
portions of the particle as the thickness of the ice 
layer increases [3]. This, combined with the facts 
that particles will be rotating, have high surface 
tension, and quickly rising internal pressure, 
allows us to assume a spherical shape for the 
particle. 


ENVIRONMENT MODEL 


The energy budget for liquid or solid water 
particles in LEO consists of four major 
components: heat flux from the Sun plus Sun flux 
re -radiated from Earth atmosphere, infra-red heat 
flux from the Earth, blackbody radiation to deep 
space, and sublimation. Heating due to radiation 
from Earth’s lower atmosphere is extremely 
small, because water vapor in the atmosphere 
absorbs pretty much in the same spectral regions 
as solid or liquid water particles in the orbit. The 
re -radiated spectrum fonn higher in the 
atmosphere has the most effect on water drops or 
ice particles in orbit. Heating from molecular 
collisions is easily shown to be negligible at 
orbital altitudes and thus is not accounted for in 
this study. Using above mentioned assumptions, 
the energy balance for the particle can be 
represented by the following formula: 


n3 dT d rop _ 

3 R d’-opP C P dr 


drop 7 B £ V drop ^ Space 


) 


where 

Qin drop ’ Sunreal + Q Earth 

7 X^Rdrop ’ A)Q Albedo ^ Q atm. heal 


( 3 ) 


( 4 ) 


LONG WAVE RADIATION ABSORPTION 

Water droplets and ice particles have very similar 
absorption characteristics with more intense 
absorption in the long wave region. It was shown 
by Kondratyev et all [4] that in the Earth infra-red 
spectral region, the main portion of which is 
localized at wavelengths from 4 to 40 p, almost all 
energy is completely absorbed by the atmosphere 
and re-radiated, leaving two transparency 
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windows in 8 to 12 p and 16 to 24 u. We have 
calculated from the energy spectrum and real and 
imaginary indices of refraction over that range 
that 7% of the flux is reflected from a crystalline 
ice particle. For our model we calculated particle 
absorptivity to be %{R drop , f) = 0.93 of Earth IR. 

For modeling purposes and in order to account for 
the spectral intensity of the scattered infra red, 
Earth is assumed to be a blackbody with a lower 
than usually assumed temperature at 25 5K. Figure 
2 shows black body emission curves for Sun and 
Earth. 


90 



Wavelength (X) in pm 

Figure 2 Black Body Emission Curves of the 
Sun and Earth 


The view factor is calculated for the altitude of the 
particle as it rapidly deorbits: 


1-Jl- 


Earth 


' REarth+S O Y 

V ^ Earth + Alt ^ 


(5) 


where 50km was chosen to be a reference 
(equilibrium) altitude for a black body 
temperature of 255 K, where IR spectrum is 
constantly absorbed and reemitted below this 
altitude. 

By calculating the Earth view factor we can 
calculate long wave radiation received from Earth 
at a specific altitude: 

Q Earth 

4 Xz(Kr 0p ’^h B R L P ( T drop ~ T L,h V Earth ’ 


SUNLIGHT ABSORPTION 

For particles larger than 100 microns in 
diameter, both Q Albedo and Q sunmerage are 

calculated using geometric optics, analyzing 
the path length of each wavelength within the 
sphere as a multiple of its absorption length. 


The total energy flux from the Sun, with a 
blackbody spectrum of 6000K, averages to 
1363 Watts/m 2 in Near Earth Orbit (NEO). 

Daily variations due to the eccentricity of 
Earth’s orbit can be calculated as: 


Q 

Q 


Sunreal 


Sunaverage 


a 


r 

0.034 cos 

V 


360/V 

365.25 


x 

/ 


(7) 


Although yearly variations of solar radiation 
are within 5% of average, they were easy to 
account for and thus were used in the model. 
These variations may play a bigger role in the 
future modeling, when the effects of 
surfactants, impurities and of surface 
roughness are considered, causing absorption 
of solar flux to be a more dominant fraction of 
the heat balance [5]. 

Below 100 micron diameter, in optics it is 
important to consider the effects of Mie 
scattering. This phenomenon causes an 
increasing portion of the light to be scattered 
as the drop radius approaches the comparable 
wavelengths. Thus, the heat input slows, and 
the sublimation rate rapidly falls off. 

Interesting Hypothesis: Droplets smaller than 
100 microns are not considered in this 
study. However there is a potential that 
due to vanishingly small sublimation rates 
and infinitesimal ballistic number, they 
may survive atmospheric impact and seed 
the highest reaches of the atmosphere as 
tiny pure ice crystals. Orbital speeds and 
atmospheric densities can decelerate a 
sub-micron particle in seconds to thermal 
speeds. This phenomenon, if verified, 
may partially explain the recent 
proliferation of noctilucent clouds, which 
were extremely rare before this century 
and which are becoming increasingly 
common in the space age. Noctilucent 
clouds appear at 80- 100km altitudes 
above Earth, and reflect sunlight long 
after dusk and long before sunrise. Such 
clouds need both concentrations of water 
vapor at high altitudes (reasonably 
common) and a large number of seed 
particles at 80- 100km altitude 
(unexplained and previously quite rare, 
occurring usually only after massive 
volcanoes). The tightly clustered de-orbit 
of billions of sub-micron particles in 


4 


spacecraft water vent operations may seed 
such clouds. 


Numerous attempts have been made to simplify 
the absorption process of sunlight in small 
droplets as a simple fraction of incident flux. 
Unfortunately, there are many nonlinear effects 
that drive this process, which must be integrated 
over all wavelengths over the entire illuminated 
surface in a range of droplet sizes. Ice exhibits a 
different real and imaginary index of refraction 
for wavelength in the solar spectrum. The 
combination of these indices determines the 
amount of reflection (which varies with 
polarization, but is always averaged 50:50 
between vertical and horizontal polarization over 
a sphere), the extent of refraction, and the fraction 
of energy in the refracted portion that is ultimately 
absorbed [6]. The absorption coefficient can be 
written in a form of the equation: 


x( R Jro P ’^)= R e- 


'^drop(K 

M m] +2) 


( 8 ) 


where 

m = n - ik (9) 

(m isa complex index of refraction). 

The real and imaginary indexes of refraction vary 
greatly by wavelength as shown in Figure 3. 



Figure 3 Real and imaginary indexes 


Our study shows that heat flux from the Sun Q stm 

absorbed by a crystalline spherical ice particle 
(with no impurities or bubbles) is negligible in 
comparison to other heat fluxes affecting the 
particle. For the studied range of particle sizes, the 
absorbed solar energy varied from 25 to 45 W/m 2 , 
which is only 1.8% to 3.3% of the incident solar 
flux, and a small fraction of the absorbed infra red 
from the Earth. The calculated values of absorbed 
fraction of solar flux versus drop size are 


presented in Figure 4, along with simplified 
empirical fit that was useful in computer 
calculations of ballistics of particles <2.5 mm 
radius. As the particle becomes smaller, the 
fraction of the sunlight absorbed gets smaller, due 
to decreasing path length for absorption of the 
different spectra. This effect is not seen in Earth 
IR absorption, due to the very short absorption 
length of IR in water relative to the smallest drop 
diameters. 


[—♦—Calculated Point —■—Empirical Fit | 



Figure 4 Sun absorption in crystalline sphere 


EARTH ALBEDO 

The term albedo (Latin for white) is commonly 
used to describe the overall average reflection 
coefficient of an object. The Earth’s albedo is 
composed of sunlight reflected from the surface of 
the planet and scattered from the top layers of the 
atmosphere. Fleat flux from reflected sunlight 
Q Albedo is absorbed by the crystalline clear 
spherical ice particle the same way as direct 
sunlight Q sunrea i . The difference between the 

absorption of Earth albedo and of direct sunlight 
is in the view factor of the sphere, and the incident 
intensity. Since the ice particle moves at very high 
velocity in the LEO and environment changes 
very fast we can use the average Earth albedo of 
0.3 throughout the illuminated part of the orbit. 

Heat flux reflected fonn the Earth can be 
calculated using the following fonnula: 

Q Albedo — 0. 3%( R dmp , A. 'jQsunreal ( 1 0) 


BALLISTIC COLLISIONS 

The energy added to the particle from direct 
elastic collision with a rarified gas can be only a 
fraction of the incident energy flux. 
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Q. 


atm. heat ^ therm 


^ R drop P ataX 


( 11 ) 


Since atmospheric densities at ISS altitudes are on 
the order of 10 11 to 10 -12 kg/m 3 and velocities 
are less than 8 km/sec, the maximum energy 
deposited to the particle cannot be greater then 6 
Watt/m 2 , even assuming that all energy was 
thermalized, which is not true. Most collisional 
energy is imparted in the form of deceleration. 
Thus for the modeling purposes Q baUistic is 
negligible compared to other energy fluxes, and is 
assumed to be equal to zero. 


ENERGY LOSS 


RADIATION TO SPACE 

The view factor to space is calculated as: 

(P Space =\~<P Earth (l 2 ) 


equilibrium temperature in minutes if external 
heat fluxes do not change. 

The heat of sublimation can be calculated using 
the following fonnula: 

jp 

Q sub =LA—, (15) 

dr 

The rate of sublimation 

—(T) = a(T )P(T (16) 

dr \2nkT 


depends on the equilibrium vapor pressure of 
liquid and solid water. Sublimation rates per unit 
exposed area are proportional to this pressure and 
temperature of the drop [8, 9], 

Equilibrium pressure for pure ice is 


P(T ) *3. 2e\0 n exp 


6110 


Pa 


(17) 


then long wave radiation emitted to space can be 
calculated as 

Qspace = 4 ^B R drop ( T drop ~ T Space }p Space > ( 1 ■ 3) 

SUBLIMATION 

In the process of formation particles will rapidly 
start loosing their energy through the process of 
evaporation and then, when they subsequently 
freeze, through the process of sublimation. This 
sublimation becomes the principal mechanism of 
energy loss. Sublimation is the process whereby a 
solid substance changes its state skipping a phase, 
going from solid to vapor. This is the only 
possible path for ice below 0.006 atm, such as the 
particle experiences in space. Figure 5 shows 
triple point chart for water [7]. 



Temperature 

Figure 5 Triple Point Chart 

The initial temperature of the drop does not affect 
the final equilibrium temperature. The drop 
reaches its freezing point in seconds and its 


PARTICLE DECAY RATE 

Once particles are formed they will start slowing 

Area 

down linearly proportional to their ratio. 

mass 

We took C D - 2, since we assumed particles to 
have spherical form. While decaying, particles 
will sublimate, constantly changing their ballistic 
number which will expedite particle decay rate. 


dF = l c 
dr 2 


A v2 
~ Pa, J 
m 


(18) 


Where — is the ratio of area of the ice particle to 
m 

the mass of the ice particle: 


A_ 

m 


7 tR 


drop 


, P drop drop 


3 1_ 

4 Pdrap R 


drop 


(19) 


If we integrate sublimation rate of the drop over 
one orbit period we will see later that we can 
approximate R d (r) to linearly decrease with 
time: 


R drop( T ) = R o ~ k P ( 2 °) 

where R a is the initial ice particle diameter, and k: 
being a very small constant sublimation constant, 
on the order of nanometers per second. 

The change in altitude with time is differentially 
linked to the change in velocity over time. Since: 
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V = ^ 

V ( R Earth + 

+ 7? £grt/l ) _ ^ dV 
[Alt + R Earth ) V 


( 21 ) 

( 22 ) 


Combining equations 18, 19, 21, and 22, we can 
get the equation: 

dAlt _ 3 C D I p A „„(p) ( 23 ) 

dr 2 p drop R drop ]] (Alt + R Earlh ) 


Where Alt,R drop and p Um are all varying with 

time. We have chosen to close the model with 
digital-numeric integration of equation (23), due 
to the variation of density and of sublimation rate 
with time and other nonlinearities. However, in 
the case of crystalline clear spherical particle it is 
possible to get a workable analytic solution for 
particle deorbit using equation ( 20 ) and assuming 
a time-invariant density profile that is purely a 
function of altitude. Then the solution can be 
closed by rearranging equation (23) to make one 
side completely dependent upon time, and the 
other side completely dependent upon altitude: 


i \(K~k{r) 

jo fu Al,(l> 

j-JfAlt + 

^ Earth )Paim( Alt )d( Alt) 

** C D Alt(0) 


(24) 


The left-hand side is a simple logarithm in time r, 
scaled by the constant sublimation rate k t and R a , 
the initial size of the drop. The initial radius of 
the particle R m and sublimation rate k, are the 
driving factors for the particle time in orbit. 
Exponentiation of both sides leads to a scaled time 
r being a complicated, characteristic function of 
Alt(t) and initial release altitude Alt(O). However, 
the right-hand side typically does not close 
analytically because of the mixture of powers and 
exponentials in altitude under the integral. It can 
be approximated to arbitrary precision by 
repetition if the techniques of integration of the 
product of functions. We can do similar tricks to 
detennine the phase angle of the particles ahead of 
the releasing spacecraft. This is essential when 
predicting the time and altitude of orbit plane 
crossings with spacecraft in lower orbits. 
Ultimately, the limitations imposed on the 
precision of our analysis by the need to assume 
overly simplistic atmosphere density and 
sublimation rate, and the ability to only 
approximate the integral, led us to develop full 
numeric integration tools to solve this problem. 


Our analysis has shown that atmospheric density 
variation has the most impact on particle 
ballistics. Atmospheric modeling is therefore most 
critical to a precise ballistic forecast especially for 
orbital objects with infinitesimally small B N . For a 
simple exponential function of density vs. altitude 
and slightly-varying reference scale height, 
depending upon which FI 0.7 number was chosen, 
the decay varies from ~17 hours to ~54 hours to 
decay from 400 to 250 km. This dependence is 
shown on Figure 6 . 



Figure 6 Atmospheric density effect on decay 
rate 


Because of the sensitivity to density, we have 
explored several different models of the 
atmosphere to process the ballistic forecast. 
Generally, all sophisticated models require precise 
inputs of location, altitude, F(10.7) (both averaged 
over three months, and daily values), geomagnetic 
potential, and beta angle. We must accept our 
largest uncertainty in the forecast from the 
atmospheric modeling, because the most accurate 
atmospheric and ballistic models often take fixed 
scalar inputs, which do not work in our case 
because of rapidly changing B N and altitude decay 
rate. We are in the process of building a simplified 
model of the atmosphere specifically for this kind 
of work as well as developing plans for fully 
integrated ballistic tool that incorporates real time 
atmospheric density and sublimation variations. 


ANALYSIS 

As a part of our water vent study we developed a 
small tool that calculates decay rate for the 
sublimating crystalline clear spherical ice particle, 
incorporating many small, nonlinear effects which 
are ultimately more important for contaminated 
water or for other effluents. We did integration 
using simple Runge-Kutta method simultaneously 

dT 

solving for three equations: temperature , 

dr 

dR , , . , dAlt , 

radius — and altitude change. 

dr d t 

We can point out three major phases in 
sublimation of deorbiting ice particle. In the initial 
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phase the temperature and the radius of the 
particle change relatively rapidly as the 
equilibrium is established. 

Our analysis begins a few seconds after the water 
is vented. We assume that particles will cool 
down to 273 K within approximately 6 to 8 
seconds after the vent and form 2mm diameter 
crystalline clear spheres. Within first 30 seconds 
of our calculation, the particle cools down to ~200 
K shrinking almost 1% from its original size. 
Figures 8 and 9 show the temperature drop rate 
and radius loss rate respectively. 



Figure 8 Particle initial temperature change rate 



Figure 9 Particle initial radius change 


The particle experiences tremendous energy loss 
during first 30 seconds primarily due to the 
sublimation process as the equilibrium 
temperature is reached. This expels approximately 
2% of particle’s original mass and volume. 

From the very first moment particle starts rapidly 
descending and looses ~10 meters of altitude in 
the first 30 seconds. 

After that, there is a long phase where the radial 
reduction rate is nearly constant with time, 
followed by a third final phase, when competing 
influences may make the rate change in the last 
stages. 


Following the establishment of thermal 
equilibrium, the rate of radial reduction of the 
particle is essentially constant. The radiated 
energy fluxes are all proportionate to the surface 
area, where the sublimation rate is proportional to 

, . dR dR 

the surface area times — . Thus, will 

d r dr 

remain nearly constant, while environment 
remains constant, in the very latest stages of 
orbital decay the variations in flux rate due to 
changing view factors will have larger effect. 
Figures 10 and 11 show changes in temperature 
and particle size. 



Seconds 


Figure 10 Particle temperature change 



in later stages, although Earth view factor rapidly 
increases, as the droplet shrinks, sunlight 
wavelengths are even less absorbed. This is due to 
the increasing effect of Mie scattering. Thus, 
competing influences try to increase the 1R flux 
rate (from Earth view factor) and to reduce total 
flux rate (from reduced absorption length, and 
Mie scattering). The rate of surface ablation 
changes proportionally to change of the 
viewfactor. When particle approaches ~320 km 
altitude the rate of orbital decay rapidly increases 
which leads to proportionate increase in 
sublimation rates and particle temperature 
increase. These dependencies are shown Figures 
12 and 13. 
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Seconds 


Figure 12 Particle sublimation rate 



The process of sublimation continues until the 
particle size becomes comparable to wavelength 
of incoming radiation, where it drastically slows 
down or may even completely stop because 
particle scatters nearly all external radiation. Our 
study has shown that in the worst case scenario of 
low atmospheric density (FI 0.7 — > 50) and cold 
environment (GMT day 173 and 0° P angle) 
particles can survive in LEO no longer then 59 
hours. 

Analysis was performed for 0° and 75° P angles as 
well as for aphelion and perihelion conditions. 
Using constant atmospheric parameters we see 
that annual solar variations change particle life in 
orbit only by 0.6% and maximum P angle 
variation (duration of eclipse) changes particle life 
by 8.1%. These variations do not play significant 
role in detennining particle life, compared to 
uncertainty in atmospheric density and in the 
range of particle sizes in the vented stream. 


almost two orders of magnitude. The change in 
atmospheric density using FI 0.7 number from 
minimal 50 to maximum 210 reduces particle life 
in orbit by 3 1 8%, and the best possible model of 
the atmosphere is therefore imperative. As 
mentioned before, no current high fidelity 
atmosphere model used for ballistic predictions 
can incorporate such rapidly changing properties 
and altitudes as our problem engenders, and we 
are in the process of building new customized 
models to address this problem. 

SURFACTANTS 

The term “surfactant” is a blend of "Surface 
active agent". Typically surfactants are organic 
compounds that contain both hydrophobic groups 
and hydrophilic groups, although almost any 
contaminant has an effect. They are typically 
soluble in both organic solvents and water. 

There is a lot of discussion in the literature on the 
effects of organic hydrocarbons on water 
evaporation and the ability of them to form 
monolayers that can significantly reduce 
evaporation. Water drop evaporation retardation 
by monolayers is possible because of their ability 
to assume more closely packed molecular 
orientation as the drop evaporates and the surface 
area contracts. Our assumption is that similar 
process will take place in space when a particle 
starts sublimating. We believe that it is possible 
for accumulating surface impurities to form a 
layer or coating on the surface of the liquid drop 
during the phase of ice particle formation and 
reduce the rate of sublimation [10, 11, 12]. The 
schematic of boundary layer formation is shown 
on Figure 14. 



Figure 14 Surface concentration 


It is clear that the smaller particle gets, the less its 
temperature deviates from equilibrium 
temperature. 

The atmospheric properties on the day of venting 
are crucial to the ballistics of these very light 
particles. Generically, the rate of decay increases 
as the particle gets smaller, altitude drops, and 
atmospheric density increases. From 400 km 
down to 250 km atmospheric density increases by 


The rate of evaporation of contaminated particle 
can be represented by the following formula: 


dR 

dr 


D 

Rdrop 


(c 0 



(19) 


As discussed above, the evaporation rate should 
remain constant if the external environment does 
not change. However the rate of evaporation 
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reduces as the concentration of surfactants at the 
outer layer of the particle increases. Figure 15 
shows the dependence of surfactants of 
evaporation rates. 



Time minuts 

Figure 15 Evaporation rate reduction 


The curvature of the drop does not play a large 
role in changing evaporation rate unless particles 
get too small, and since we are not concerned 
about very small particles we can neglect this 
effect: 


In 


c 


o 


2 oM 

KopRTp 


( 20 ) 


The rate of evaporation of the particle containing 
surface impurities can be represented as: 


dR dro P _ D (c 0 ~c g ) 
dr R drop u KR 2 0 ’ 

D 3 


( 21 ) 


where K is found empirically. 
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Figure 16 Contamination size distribution 

The presence of humans on board detennines 
waste water composition and existence of 
different types of hydrocarbons including long 
chain alcohols, which act as effective surfactants. 
Figure 17 shows ISS waste water composition 
[13]. 
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Figure 17 Waste water composition 


Water samples were collected through several ISS 
expeditions and analyzed by US and Russian 
specialists. 


It is still uncertain how fast the outer layer is 
fonned and if it is carried away by initial process 
of particle fonnation. If it is formed within the 
time it takes for the particle to freeze, it is possible 
that the monolayer might be expelled in a kind of 
purifying step, in the violent evaporation phase. If 
not, it may remain as an inhibiting seal on the 
surface of the particle. 

WASTE WATER CONTAMINATION 

Currently there are two manned space vehicles 
designed to release waste water into space. The 
International Space Station (ISS) and Space 
Shuttle are each capable of releasing ~60 liters of 
waste water at one vent. The size of the cloud and 
the process of particle formation are similar for 
both vehicles. The major difference between the 
two is the duration of particle survival in LEO due 
to different impurities contained in vented water. 
We have conducted a study of ISS waste water 
and determined that it contains mechanical 
impurities that differ in size from 25 to 400 
microns. Figure 16 shows size distribution of 
impurities. 


Space Shuttle waste water will contain more 
organics in it than ISS condensate based on the 
design of the system. We believe that this will 
significantly slow sublimation of Space Shuttle 
waste water particles. At this moment we cannot 
definitely say by how much the process of 
sublimation will be retarded. 

CONCLUSIONS 

As a result of this study we were able to show that 
particle survival in LEO primarily depends on the 
density of the atmosphere at the time of the vent. 
In most cases of average or high atmosphere 
particles will reach safe altitudes before they 
sublimate to safe sizes. Uncertainties in 
atmospheric density, bubble inclusion, and the 
presence of surfactants all add to a large 
uncertainty band when forecasting exact 
trajectories for a family of candidate drop sizes in 
the vented plume. In the worst case scenario of 
low atmosphere and cold environment, the largest 
particles will sublimate to sub-millimeter sizes 
within 2.5 days, thus presenting a potential threat 
to lower-orbiting spacecraft. Based on the results 
of this study and to mitigate the risk of ice 
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particles from ISS waste water damaging Shuttle 
tiles we developed a flight rule that restricts any 
waste water vents by ISS three days prior to any 
Shuttle launch. We are in the process of refining 
the analytic tools to narrow the predicted range of 
altitudes and phase angles of the vented stream of 
particles, such that planned ventings can be timed 
to preclude intersection of the vented particles 
with any high-value spacecraft in crossing orbits 
below the vent. A new class of ballistic 
forecasting tool must be developed to accomplish 
this, accepting dramatically carrying ballistic 
number and rapid altitude decay in a precise 
density profile that varies with location, altitude, 
and daily solar and geomagnetic parameters. 

The effect of surfactants on sublimation of ice 
particles should be studied further to evaluate the 
effect on sublimation rates and reduce the 
potential hazard of catastrophic event to a lower 
flying vehicle. 
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Problem statement 


• In vacuum large diameter streams with high vapor 
pressures have a tendency to burst, forming an 
uncontrolled cloud of particles which subsequently 
freeze and eventually sublimate 

• The behavior of the stream and the formation of 
the particles in vacuum depend on several factors 
such as vapor pressure, viscosity of the vented 
fluid, dissolved gas content, surface tension, 
surfactants and the density 

• In this study we attempted to capture the behavior 
of the particle in the Low Earth Orbit, its 

ht sublimation and deorbit process 
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Assumptions 

• Particle Model 

■ Spherical crystalline clear ice particle 

■ No impurities or organic solvents. 

■ Clear surface with no roughening 

■ Homogeneous with no dissolved bubbles 

■ Background temperature for Earth 255K at 50 km reference altitude 

■ Background temperature of Space 3K. 

■ Heat flux absorption by the ice particle was calculated from complex 
indices of refraction 

> Earth IR = 0.93 

> Direct Sun flux = 0.27 to 0.06 

> Earth albedo = 0.3 


Atmospheric model 

■ Particles remain in circular or nearly circular orbit 

■ Initial venting altitude = 400 km. 

The coefficient of drag = 2, characteristic of very smooth spherical 
irticles. 
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Short wave radiation absorption 
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Black Body Emission Curves of the Sun and Earth 
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Conclusions 


Particle survival in LEO primarily depends on the density of the atmosphere at the 
time of the vent 

Uncertainties in atmospheric density, bubble inclusion, and the presence of 
surfactants all add to a large uncertainty band when forecasting exact trajectories 
for a family of candidate drop sizes in the vented plume. 

Based on the results of this study we developed a flight rule that restricts any 
waste water vents by ISS three days prior to any Shuttle launch 

We are in the process of refining the analytic tools to narrow the predicted range of 
altitudes and phase angles of the vented stream of particles to preclude 
intersection of the vented particles with any spacecraft in crossing orbits below the 
vent 

The effect of surfactants on sublimation of ice particles should be studied further 
to evaluate the effect on sublimation rates and reduce the potential hazard of 
catastrophic event to a lower flying vehicle 


